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The electrostatic potential above the Abrikosov vortex lattice, discussed earlier by Blatter et al. 
[PRL 77, 566 (1996)], is evaluated within the Ginzburg-Landau theory. Unlike previous studies 
we include the surface dipole. Close to the critical temperature, the surface dipole reduces the 
electrostatic potential to values below sensitivity of recent sensors. At low temperatures the surface 
dipole is less effective and the electrostatic potential remains observable as predicted earlier. 



I. INTRODUCTION 

A boundary between the normal and superconducting 
states is characterized by two length scales. First, the 
local fraction of superconducting electrons is measured 
by the Ginzburg-Landau (GL) wave function ip which 
changes on the scale of the GL coherence length £. Sec- 
ond, there is a magnetic held B screened on the scale of 
the London penetration depth A. This picture is common 
to type-I superconductors, in which the superconducting 
state is nearly separated from the normal state, and to 
type-II superconductors, in which the normal state is dis- 
persed into individual lines called vortices. 

In extreme type-II superconductors, where £ <C A, the 
spatial shape of the superconducting fraction provides a 
sharper image of vortices than the spatial dependence of 
the magnetic field. Blatter et al. 1 proposed to observe 
the space modulation of the superconducting fraction via 
the electrostatic field that is expected to develop above 
a surface of superconductors. Their estimate of the elec- 
tric field created by the Abrikosov vortex lattice predicts 
values well observable by recent experimental tools. 

The idea of Blatter et al. is as follows. According to 
theoretical predictions, 2-4 the space modulation of the 
superconducting gap A induces a charge transfer so that 
an electrostatic potential (called the Bernoulli potential) 
develops inside the superconductor. Since the Bernoulli 
potential (f> is a function of the square of the gap, <fi can 
be used to observe | A| 2 . The GL wave function is linearly 
proportional to the superconducting gap, i.e., </> can be 
used to observe \^\ 2 . 

The Bernoulli potential cannot be detected inside the 
superconductor, it leaks out from the surface, however. 
Its detection outside might be possible by scanning force 
microscopy or using the Kelvin capacitive pickup with a 
single-electron transistor as a sensor. 1 

As far as we know, such an experiment has not been 
performed yet, but it is under preparation. It is likely 



that the experimental setup will first be tested on con- 
ventional superconductors. 5 In this paper we show that 
close to the critical temperature the electrostatic poten- 
tial above the surface is strongly reduced by the surface 
dipole, which arises due to unbalanced pairing forces. 6 At 
lower temperatures the surface dipole is not so effective 
and the electrostatic field reaches observable amplitudes. 

The paper is organized as follows. In the next sec- 
tion we specify the assumed experimental situation and 
the physical picture of charge transfer contributing to 
the expected electrostatic potential. In Sec. Ill we de- 
rive a relation between the GL wave function and the 
electrostatic potential and provide a simple estimate of 
the potential for the magnetic field far from the upper 
critical field, i.e., for the limit of separated vortices. In 
Sec. IV we discuss numerical results, and in Sec. V we 
summarize. 



II. BASIC ASSUMPTIONS 

Let us first describe the experimental situation we as- 
sume in our discussion. A superconducting film fills the 
layer — L < z < 0. This film is thin on the characteristic 
scales of the GL theory, L <C £, A, but it is thick with 
respect to the BCS coherence length, L ^> £ - The mag- 
netic field B||z penetrates the superconductor in form of 
vortices. The electrostatic potential will be scanned close 
to the surface at z scan > 0. 

Now we turn to the underlying physical picture. The 
electrostatic potential leaking out of the superconductor 
is generated by charges that can be sorted into three 
groups: A. the bulk charge, B. the surface dipole, and 
C. the surface charge. These contributions are discussed 
in individual subsections. 
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A. Bulk charge 

The bulk charge covers a transfer of electrons from the 
inner to the outer regions of the vortices. There are var- 
ious forces taking part in this transfer. First, electrons 
rotate around the vortex center so that the inertial force 
acts in the centrifugal direction. Second, the magnetic 
field pushes electrons via the Lorentz force - also pointing 
outward. Third, the energy of Cooper pairs is lower than 
the energy of free electrons, therefore unpaired electrons 
in the vortex core are attracted towards the condensate 
around the core. The resulting force again points out- 
ward. These forces deplete the electron density in the 
vortex core, creating the Coulomb force which balances 
all the other forces. 

Due to these various contributing mechanisms, there 
is no single characteristic length scale of the charge mod- 
ulation. The shortest scale is the GL coherence length 
£ reflecting the pairing forces. The contribution of the 
inertial and the Lorentz forces change on the scale of the 
London penetration depth A. The long-range periodicity 
is of course enforced by the Abrikosov vortex lattice. 

The bulk charge has been evaluated within various ap- 
proximations mostly covering only some of the acting 
forces. Studies based exclusively on the Lorentz and the 
inertial forces have been performed within the classical 
picture of the superconducting fluid 7 and later with its 
quantum form. 8 

In the last decade, the charge transfer in vortices has 
been derived from the electron-hole non-symmetry of the 
density of states at the Fermi level. 3 ' 4 Later it was rec- 
ognized that the electron-hole non-symmetry effects are 
identical to forces due to the pairing energy. 9 ' 10 In result 
there are two scientific dialects for the pairing forces. We 
prefer the original one. 

Approaches that cover all the above listed forces are ei- 
ther phcnomenological or microscopic. The phenomeno- 
logical studies 11 use the theory of the GL type. The 
microscopic studies 12-16 are based on the Bogoliubov- 
de Gennes theory. It is worth to mention mention that 
in the vicinity of T c the electrostatic potential obtained 
from the microscopic studies agrees with the result of the 
GL theory. 12 " 16 

The microscopic theory is superior to the others as it 
covers all important contributing mechanism in a unified 
way 16 and its region of applicability is not restricted to 
the vicinity of T c or to small gradients of the supercon- 
ducting gap. On the other hand, such complete treat- 
ment is complicated and numerically demanding. Here 
we use the simpler GL theory. 

B. Surface dipole 

The second kind of charges determining the electro- 
static potential is a surface dipole due to which the po- 
tential has a finite step. A spatially resolved profile of this 



step is not known. It is expected that the dipole is located 
near the surface inside the superconductor. In analogy 
with the space profile of the superconducting gap, one 
can speculate that the width of the dipole is similar to 
the BCS coherence length £ - The present treatment is 
limited to temperatures close to T c , where £o <C £■ Ac- 
cordingly, we assume the width of the surface dipole to 
be infinitesimal. 

So far, the role of the surface dipole in superconductors 
is not fully clarified. There is an important experimen- 
tal experience with the electrostatic potential above the 
surface in the Meissner state, however. Precise measure- 
ments of the electrostatic potential made by Morris and 
Brown 17 with the help of the capacitive Kelvin method 
have shown a surprising result - all contributions of the 
pairing forces 9 ' 10 are canceled by the surface dipole. The 
experiment of Morris and Brown thus indicates that the 
surface dipole has to be taken into account. 

Based on the assumption that the surface dipole results 
from the surface depression of the BCS gap, we have de- 
rived in Ref. 6 the local value of the dipole as a function 
of the GL wave function at the surface. The internal 
electrostatic potential and the potential step due to the 
surface dipole add. The resulting observable surface po- 
tential 

e^o = -- (1) 
n 

is given by the free energy per electron. 6 Here n is the 
total density of pairable electrons (the total density) and 
/el is the density of the electronic part of the free energy. 

The surface potential (1) follows from general thermo- 
dynamic assumptions and it can be implemented within 
different approximations of the free energy, compare 
Ref. 6 with Ref. 18. In our treatment the free energy 
is evaluated from the GL theory. 

C. Surface charge 

Finally, there is the surface charge distributed on the 
scale of the Thomas- Fermi screening length, Atf, from 
the surface. Since Atf is much shorter than the GL co- 
herence length £, the BCS coherence length £ , and the 
London penetration depth A, we use the limit Atf — * 
and treat the surface charge as an ideal two-dimensional 
surface charge. The surface charge represents the utmost 
layer of the superconductor. 

The surface charge simplifies the construction of 
the electrostatic potential inside and outside of the 
superconductor. 1 Neglecting contributions of the order 
of A TF /^ 2 , one can simply evaluate the potential at the 
inner side of the surface, </>(0_), and match it with the 
potential outside, <j>(0+) = 0(O_). The potential outside 
decays far from the surface to its mean value, <p(z) — > ((f)) 
for z — > oo. Due to this asymptotic condition, the po- 
tential is fully specified by its value at the surface. In 
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our case the potential at the surface includes the sur- 
face dipole, 0(0+) = 0(O_ ) = <po, where 4>o is given by 
formula (1). 

The matching of the inner and outer potentials is sim- 
ple in the two-dimensional Fourier representation, 

0(K) = ^jf drfo(r)cos(Kr), (2) 

where r = (x, y) and K are discrete momenta that have 
to be selected according to the structure of the Abrikosov 
vortex lattice. The area of the elementary cell is given 
by the mean magnetic field B and the elementary flux, 
<f>o = Bfl. One obtains the potential at any distance 
z > from the surface, 

0(r, z) = (cj>) + ^ 0(K) c-l K l* cos (Kr) , (3) 

from the Fourier components (2) and the mean value 
(4>) = / dr o (r). It is easy to check that the potential 
(3) satisfies the Poisson equation. 

The potential or its gradient at a finite distance will 
naturally be necessary for interpretations of future mea- 
surements. On the other hand, we feel that this tech- 
nical step does not bring any new insight to the prob- 
lem. In our discussion we will focus on the surface value 
0o(r)=0(r,O). 

III. SURFACE POTENTIAL WITHIN THE 
GINZBURG-LANDAU THEORY 

With respect to electrostatic phenomena it is advan- 
tageous to introduce the GL theory in the formulation 
proposed by Bardeen. 19 ' 20 The free energy 

/el = \lT 2 + (-iftV - e*Af i, 

2 2m 

n 2 V n 

is composed of three terms: the free energy of Gorter and 
Casimir (last two terms), the kinetic energy in the quan- 
tum form (second term), and the subtracted free energy 
of the normal state (first term). Here, 7 is the linear 
coefficient of the specific heat. The condensation energy 
determines the critical temperature as £ con = jlT 2 . 

For temperatures close to T c the superconducting frac- 
tion is small, 2 ^ <C 1. Expanding the square root in 
this small value one arrives at the GL free energy 

/el = HW - e*A) 2 4, + a\ip\ 2 + i/3H 4 . (5) 

Note that the kinetic energy is not in the form proposed 
by Ginzburg and Landau 21 . In particular, it becomes 
negative at the vortex core. The volume integral over 
this kinetic energy differs from the GL form by a surface 



contribution which is identically zero if the GL boundary 
condition is satisfied, see the discussion in Rcf. 22. 

From the limit ► one obtains the GL param- 
eters a = -^7 (T 2 - T 2 ) and (3 = ^jT 2 . Below we 
compare the surface potential with the bulk potential. 
The bulk potential includes pairing forces proportional 
to density derivatives ^ and The temperature de- 
pendence is essential in this point as 7^ 0. After we 
have evaluated derivatives, we can use the usual limiting 
parameters of the GL theory, a = ——^T c (T c — T) and 

f> = 

To complete the description of the GL free energy we 
recall that the charge of the Cooper pair equals twice 
the electron charge, e* = 2e. The mass m* — 2m de- 
pends on impurities and has to be fitted, e.g., from the 
GL coherence length using £ 2 to* 7 (T 2 - T 2 ) = nh 2 . For 
Niobium we use T c = 9.5 K and 7 = 719 J/(m 3 K 2 ) 
giving e con = 1.6 10 4 J/m 3 . The electron density is 
n = 2.2 1 28 /m 3 , so that the condensation energy per 
particle is e con /n = 4.59 ^cV. For pure Niobium the ef- 
fective mass is m = 1.2 m e and the GL parameter is 
k = 0.78. 



A. The surface potential 

To proceed we have to derive the GL equation for the 
GL wave function For the assumed thin layer, induced 
currents are proportional to the layer width L and they 
are negligible in the limit L <C A. The vector potential 
A thus has the same value as in the absence of the su- 
perconducting layer and B is practically constant in thin 
layers 23 . 

From the minimum of the free energy we arrive at the 
GL equation 

(-iTN - e*A) 2 2 

_ 1[> = -ai>- PWi). (6) 

This equation is solved numerically with the help of the 
iteration procedure described in Ref. 24. 

Before we enter a discussion of numerical results pre- 
sented in the next section, it is useful to express the sur- 
face potential as a function of the amplitude of the GL 
wave function. To this end we multiply (6) by 4> which 
yields 2^4>{-iKV - e*A) 2 ?/> = -a|^| 2 - /3|V>| 4 - If we 
substitute this relation into the free energy (5), we find 
that the free energy is the negative value of the quar- 
tic term j c \ = — i/3|i/'| 4 . The surface potential (1) thus 
attains the simple form 

e^o = ^/3M 4 . (7) 

At this point we can stress the first effect of the surface 
dipole. According to (7), the electrostatic potential is 
proportional to |V'| 4 . In contrast, the internal potential 
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used by Blatter et al. 



e0Bi 



n On 



(8) 



is proportional to \ip\ 2 . The different order of the ampli- 
tude does not ruin the basic idea of the proposed exper- 
iment - the electrostatic potential at the surface can be 
used to monitor the GL wave function. 



B. Upper estimate of the amplitude 

It is possible to establish an upper estimate of the am- 
plitude of the electrostatic potential. The minimum of 
the potential (7) is at a vortex center, where ecfio — » 
since \ip\ 2 — > 0. 

The potential e<f>Q reaches its maximum somewhere be- 
tween vortices. The magnetic field in the Abrikosov vor- 
tex lattice suppresses the amplitude of the GL wave func- 
tion compared to its value xp 2 ^ = — jj in the non-magnetic 

state. Therefore from \ip\ 2 < Voo one obtains an upper 
estimate of the maximum as 



_L^1 

2n 



-con 

n 



4(l-t) 2 



(9) 



Here we have used a = (— 4e con /n)(l — t), where t = 
T/T c , to highlight that close to the critical temperature 
the amplitude of the potential vanishes as (1 — t) 2 . 

To illustrate how the surface dipole changes the sur- 
face potential we compare the full potential (7) with the 
internal potential (8). For the amplitude of the internal 
potential (8) one has the upper estimate 



< - 



7T C dT c a 
n dn (3 



£con91nTc -8(l-t). (10) 



n 



d Inn 



For Niobium one finds from the McMillan formula the 
value \"yT c ^- = 6.78 /tteV, see Ref. 11. This corresponds 



to 



dlnT c 
dlnn 



= 0.74. Numerical factors of the surface poten- 
tial (9) and the internal potential (10) are thus quite 
comparable. What makes the difference is the tempera- 
ture dependence. 



C. Lower temperatures 

The above upper estimate of the potential amplitude 
is based on the zero free energy at the vortex center and 
the finite free energy in the non-magnetic state. These 
values can be easily estimated also at temperatures far 
from T c . Indeed, as long as the free energy is adjusted to 
vanish in the normal state, it has to be zero at the vortex 
center, where the superconducting condensate vanishes. 
The free energy in the non-magnetic state can be taken 
from the thermodynamic measurements in the form / e i = 



— £ C on(l — t 2 ) 2 . Accordingly, the upper estimate of the 
amplitude of the potential is 



e0o<— (l-t 2 ) 2 . 
n 



(11) 



This estimate applies to any temperature. Of course, it 
approaches (9) for t — > 1. 

Similarly we can use the relation between the GL wave 
function and the superconducting density 2 \ip\ 2 = n s to- 
gether with the phenomenological law n s = n (1 — t 4 ) to 
obtain the estimate of the internal potential 



dlnT r 



d Inn 



2 (l-t 4 ). 



(12) 



For t — ► 1 this estimate goes to (10). 

One can see that for low temperatures, t < 1, esti- 
mates (11) and (12) yield comparable amplitudes of the 
potential. This indicates that at low temperatures the 
surface dipole is less effective in reducing the surface po- 
tential. 



IV. NUMERICAL RESULTS FOR NIOBIUM 

Now we discuss the electrostatic potential with the 
help of numerical results. To be specific, we use material 
parameters of Niobium. 

Figure 1 displays the GL wave function u) = \ip\ 2 /tp"^, 
where Voo is the GL wave function in the absence of 
the magnetic field, 2-0 2 < , = n(l — t 4 ). The profile of u> 
is compared with the surface electrostatic potential 4>o 
according to equation (7) near the critical temperature, 
T = 0.95T C . 

In the first row we show the result for Niobium with 
k = 1.5. The low magnetic field B — 0.06 B c2 already 



falls into the limit of isolated vortices, because vortex 
cores are well separated and the superconducting fraction 
between them reaches its non-magnetic value with less 
than 1% difference. The amplitude of the electrostatic 
field is thus well described by the estimate (11). 

Compared to the superconducting fraction, the elec- 
trostatic potential is much flatter in the center of the 
vortex. This feature directly follows from formula (7). 
At the vortex center, x 2 + y 2 = r 2 — > 0, the supercon- 
ducting fraction quadratically vanishes with r, oj cx r 2 . 
According to (7) the potential vanishes there with the 
fourth order of r, O oc r 4 . We believe that this will be 
possible to be observed in future measurements. 

The second row in figure 1 with k = 0.78 corresponds 
to pure Niobium. We choose the magnetic field close 
to the upper critical field B = 0.7818 B c2 , when the 
superconducting fraction is suppressed to less than 1/3 
of its non-magnetic value. In this regime the GL wave 
function is well approximated by the asymptotic solution 
due to Abrikosov. 25 
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B/B c2 (T)=0.06 CO [nV] 




FIG. 1. The relative superconducting fraction uj = \ip\ 2 /ip^ (left column) and the electrostatic potential (f> (right column) 
at the surface of the superconductor with the Abrikosov vortex lattice. The distances are normalized to the vortex separation. 
The temperature T = 0.95 T c is used for both cases, the GL parameter n and the magnetic field are specified for each row. 
Material parameters are of Niobium. 
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A. Estimate from the Abrikosov solution 

For B close to B c2 , the amplitude of the GL wave 
function undergoes rapid changes while its shape remains 
nearly unchanged. As in Fig. 1 we express the supercon- 
ducting fraction, 2\tp\ 2 /n, in terms of its non-magnetic 
value, 2\ip O0 \ 2 /n — 1 — t 4 as the relative superconducting 
fraction w = |-!/>| 2 /|V'oo| 2 - 

For the fixed shape of the GL wave function, the mean 
(ui) = J dr uj and its fluctuation (uo 2 ) = J dr uo 2 
are proportional as {ui 2 ) = /3a (w) 2 , where the Abrikosov 
coefficient for the triangular lattice is /3a = 1.16. Taking 
1 — b with b = B/B c2 as a small perturbation, one finds 
that the free energy (5) has the minimum when (uj 2 ) = 
(w) (1 — b). The mean and the fluctuations thus read* 
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PA 



<- 2 > = ^ (1 
PA 



(13) 



With the help of the mean values (13) we can express 
the mean values of the electrostatic potentials. The sur- 
face potential (7) with the surface dipole included has 
the mean value 

(e^) = £ —(l-t 2 f^(l-bf. (14) 
n pa 

If we extend the estimate by Blatter et al. to magnetic 
fields close to the critical field, we find that it has the 
mean value 

Comparing (14) with (15) one can see that due to the 
surface dipole, the electrostatic potential vanishes close 
to the upper critical field as (1 — b) 2 , rather than 1 — 6 
without the dipole. We expect that this dependence on 
the magnetic field might be one of experimentally well 
accessible tests of the presence of the surface dipole. 



B. Effect of the surface dipole 

As already mentioned the surface dipole is responsible 
for the flat region of the potential at the center of the 
vortex. In this section we discuss the role of the surface 
dipole in more detail. 

In Fig. 2 we compare the present theory yielding for- 
mula (7) with three approximations. First, if one neglects 
the surface dipole, the surface potential becomes equal to 
the internal potential 11 



*For details of the derivation see Ref. 25. In the limit of the 
thin layer the term (uj 2 )/2k 2 disappears as it follows from the 
induced magnetic field. 
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i\> (—iKV — e*A) V 
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+ 



2n 2 



(16) 



The material parameters for Niobium d6 Q°' 1 =8.71 /xeV 
and 5§JT 2 = 3.85 ueV, arc derived in Ref. 11. 

Second, the surface dipole results from the pairing 
forces. As long as one does not account for the surface 
dipole, perhaps it is better to neglect also other pairing 
forces. In this approximation the surface potential equals 
the first term of equation (16) which covers the inertial 
and Lorentz forces. 

Third, following Blatter et al. we take a single vortex 
located at x,y = 0. Its GL wave function is modeled by 
the Clem ansatz 26 wrj 1 — £ 2 /(x 2 + y 2 +£ 2 ). The vortex 
radius £ v found fro m the minimum of the fre e energy is 
given by £ v = £y/2y/l - K 2 (^/\)/K 2 (^/\), where K 
and K\ are modified Bcsscl functions. According to the 
approximation of Khomskii and Kusmartsev 3 adopted by 
Blatter et al., we take only the second term of equation 
(16) W ith %^ « I 7 T C ^ = 6.78 ^eV. 




0.4 0.6 

x,y 

FIG. 2. Various approximations of the electrostatic poten- 
tial at the surface. The parameters T = 0.95 T c , k = 1.5 
and B — 0.7818 B C 2 correspond to the upper right plot in 
Fig. 1. The symmetric lines (solid) are cuts along the x axis, 
the non-symmetric lines (dashed) are cuts along the y axis in 
Fig. 1. The line 'present theory' is given by eq. (7), the line 
'no surface dipole' by eq. (16), the line 'inertial + Lorentz' 
corresponds to the first term of eq. (16), and the line 'Blatter 
et al.' is according to Ref. 1. 

The potentials plotted in Fig. 2 can be sorted into two 
groups. The internal potential (16) and Blatter's result 
are very similar except for some minor differences follow- 
ing from the Clem model and the neglect of The 
potential (7) with the surface dipole included, and the 
approximation by inertial and Lorentz forces are much 
smaller than the potentials from the first group. 
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-50 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

0.2 0.4 0.6 0.8 1 

FIG. 3. Detail of Fig. 2. 

To compare our potential (7) with the approximation 
by inertial and Lorcntz forces, we have expanded the 
scale in Fig. 3. The numerical result clearly shows that 
the approximation by inertial and Lorentz forces has a 
very different profile. Briefly, above the Abrikosov vor- 
tex lattice the surface dipole cancels the major part of 
the contribution of pairing forces. On the other hand, it 
is not possible to avoid the full calculation and to replace 
the action of the surface dipole by simply omitting the 
pairing forces. 

Figure 4 shows the different potentials for k = 0.78 cor- 
responding to pure Niobium of the lower row in Fig. 1. 
Most of the features are identical to the situation of larger 
k in the upper row in Fig 1 and in Fig. 2. But, in the 
internal potential (16) the GL wave function is now sup- 
pressed by the magnetic field compared to Blatter's ap- 
proach. Indeed, the Clem approximation is derived for 
the limit of low magnetic fields and thus it does not cover 
the suppression. 



^ -0.5 



-1.5 
-2 



inertial + Lorentz 




x ^\ present theory 




\ 

\ no surface dipole 




\ 
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0.2 0.4 0.6 0.8 1 



x, y 

FIG. 4. The same as Fig. 2 but with parameters of the 
lower row in Fig. 1, i.e., T = 0.95 T c , k = 0.78 and 
B = 0.78 B c2 . 
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FIG. 5. Detail of Fig. 4. 

The detail presented in Fig. 5 demonstrates that the 
full theory and the approximation by inertial and Lorentz 
forces result in very different profiles of the electrostatic 
potential. One can see that the neglect of pairing forces 
leads to the potential of much larger amplitude and the 
wrong sign. 

V. SUMMARY 

We have evaluated the electrostatic potential above 
the surface of a thin superconducting layer with the 
Abrikosov vortex lattice. It has been shown that the 
surface dipole strongly modifies the magnitude of this 
potential, in particular when the GL wave function has 
a small magnitude. This is due to the relation 0o oc \ip\ 4 , 
see (7), while without the dipole one finds 4>b\ oc l^l 2 - 

According to various mechanisms to suppress the GL 
wave function, we can outline possible cases for which the 
presented theory can be tested. At the vortex core \ip\ 2 oc 
r 2 so that <fio oc r 4 while <pB\ oc r 2 . At temperatures close 
the critical temperature, t — > 1, \ip\ 2 oc 1 — t, therefore 
4>o oc (1 — t) 2 while </>bi oc 1 — t. Finally, at magnetic fields 
close to the upper critical field, b — > 1, |-0| 2 oc 1 — 6 so 
that 0o oc (1 — b) 2 while </>bi oc 1 — b. Of course, there is 
a rich area of experimental realizations with mesoscopic 
superconducting devices. 

We would like to mention limitations of our approach. 
First, the formula for the surface potential has been de- 
rived only for the magnetic field parallel to the surface. 
According to its interpretation in terms of the pairing en- 
ergy we believe that it also applies for the perpendicular 
field, nevertheless, its validity should be tested. Second, 
the local approximation of the surface dipole requires the 
BCS coherence length £o to be much smaller than the GL 
coherence length £. This is satisfied at temperatures close 
the critical temperature, while one can expect sharper 
spatial profiles at lower temperatures. For this region of 
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lower temperature, however, our results have to be taken 
only qualitatively. 
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